EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




CERN-PH-EP/2011-034 
2013/01/19 



CMS-BPH-10-005 



Measurement of the B° Production Cross Section in pp 

Collisions at \/s = 7TeV 



The CMS Collaboratioif] 



Abstract 

Measurements of the differential production cross sections da I dpj and da/ dy B for 
B° mesons produced in pp collisions at = 7TeV are presented. The dataset used 
was collected by the CMS experiment at the LHC and corresponds to an integrated 
luminosity of 40 pb -1 . The production cross section is measured from B° meson de- 
cays reconstructed in the exclusive final state J/t/?Kg, with the subsequent decays 
]/xp — > u + u~ and Kg — > tc + tc~ . The total cross section for pj > 5GeV and \y B \ < 2.2 
is measured to be 33.2 ± 2.5 ± 3.5 yb, where the first uncertainty is statistical and the 
second is systematic. 
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Cross sections for heavy quark production in hard scattering interactions have been studied at 
pp colliders at center-of-mass energies from 630 GeV [1J to 1.96 TeV 0-31 and in p— nucleus 
collisons with beam energies from 800 to 920 GeV 0. The expected cross sections can be cal- 
culated in perturbative Quantum Chromodynamics. The comparison between data and pre- 
dictions provides a critical test of next-to-leading order (NLO) calculations [6|. Considerable 
progress has been achieved in understanding heavy quark production at Tevatron energies, 
largely resolving earlier discrepancies [7], but substantial theoretical uncertainties remain due 
to the dependence on the renormalization and factorization scales. Measurements of b-hadron 
production at 7 TeV provided by the Large Hadron Collider (LHC) EHE] represent a test at a 
new center-of-mass energy of theoretical approaches that aim to describe heavy flavor produc- 
tion mum. 

This Letter presents the first measurement of the B° cross section in pp collisions at \fs = 7 TeV. 
Events with B° mesons reconstructed from their decays to the final state J/ipKg, with ]/ip — > 
and Kg — > n + n~ , are used to measure da I dpj, da I dy B , and the integrated cross section 
for transverse momentum pj > 5 GeV and rapidity |i/ B | < 2.2, where y is defined as \ In fj^p 
E is the particle energy and pi is the particle momentum along the counterclockwise beam 

direction. As the B° and B° are indistinguishable in this analysis, both mesons are referred to 
as B° for the purposes of reconstruction and the final results are divided by two to obtain an 
average. 

The data sample collected by the Compact Muon Solenoid (CMS) detector at the LHC corre- 
sponds to an integrated luminosity of 39.6 ± 1.6 pb -1 and represents the entire 2010 dataset. A 
detailed description of the detector may be found elsewhere IIT31 . The main detector compo- 
nents used in this analysis are the silicon tracker and the muon systems. 

The silicon tracker measures charged particles within the pseudorapidity range \t]\ < 2.5, 
where f] = — In tan ( | ) and 8 is the polar angle of the track relative to the counterclockwise 
beam direction. It consists of 1440 silicon pixel and 15 148 silicon strip detector modules and 
is located in the 3.8 T field of the superconducting solenoid. It provides an impact parameter 
resolution of ~ 15 }im and a pi resolution of about 1.5% for particles with transverse momenta 
up to 100 GeV. Muons are measured in the pseudorapidity range \rj\ < 2.4, with detection 
planes made using three technologies: drift tubes, cathode strip chambers, and resistive plate 
chambers. 

Events are selected by a trigger requiring two muons without any explicit requirement on the 
muon momentum. The muon candidates are fully reconstructed offline, combining informa- 
tion from the silicon tracker and muon detectors, and are required to be within the following 
kinematic acceptance region: pl^ > 3.3 GeV for \rji'\ < 1.3; total momentum p^ > 2.9 GeV for 
1.3 < | | < 2.2; and p\ > 0.8 GeV for 2.2 < \rf^\ < 2.4. Opposite-sign muon pairs are fit 
to a common vertex to form J/tp candidates, which are required to be within 150 MeV of the 
world-average ]/tp mass Iil4|. 

The Kg candidates are formed by fitting oppositely charged tracks reconstructed with the CMS 
tracking algorithm |Tl5l to a common vertex. Each track is required to have at least 6 hits in 
the silicon tracker, a normalized x 2 < 5, and a transverse impact parameter with respect to 
the luminous region greater than 0.5 times its uncertainty. The reconstructed Kg decay vertex 
must have a normalized x 2 < 7 and a transverse separation from the luminous region at least 5 
times larger than the uncertainty on the separation. The n + n~ invariant mass m K o is required 
to satisfy 478 < m K o < 518 MeV, and the reconstructed mass distribution is found to be in good 
agreement with the world-average value l|T4| . 
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The B candidates are formed by combining a J/tp candidate with a Kg candidate. A kinematic 
fit is performed with the two muons and the Kg candidate, in which the invariant masses of 
the ]/tp and Kg candidates are constrained to their world-average values |14| . The B° vertex 
fit confidence level is required to be greater than 1% and the reconstructed B° mass must 
satisfy 4.9 < < 5.7 GeV. When more than one candidate in a single event passes all the 
selection criteria, only the candidate with the highest B° vertex fit confidence level is retained, 
which results in the correct choice 99% of the time in simulated events containing a true signal 
candidate. A total of 23 174 B° candidates pass all selection criteria. 

The efficiency of the B° reconstruction is computed with a combination of techniques using 
the data and large samples of fully simulated signal events generated by PYTHIA 6.422 [fT6]|, 
decayed by evtgen [17], and simulated by GEANT4 ||T8| . The trigger and muon-reconstruction 
efficiencies are obtained from a large sample of inclusive ]/xp — > decays in data using a 

technique similar to that described in Ref. HH, where one muon is identified with stringent 
quality requirements, and the second muon is identified using information either exclusively 
from the tracker (to measure the trigger and muon-identification efficiencies), or from the muon 
system (to measure the silicon tracking efficiency). Since the dimuon efficiencies are calculated 
as the product of the measured single muon efficiencies, a correction (1-6%), obtained from the 
simulation, is applied to take into account efficiency correlations between the two muons. The 
probabilities for the muons to lie within the kinematic acceptance region and for the B° and Kg 
candidates to pass the selection requirements are determined from the simulated events. To 
minimize the effect of the PYTHIA modeling of the p B and |y B | distributions on the efficiency 
calculation, the simulated events are reweighted to match the kinematic distributions observed 
in the data. The efficiencies for hadron-track reconstruction |20| , Kg reconstruction [21J, and for 
fulfilling the vertex quality requirement are found to be consistent between data and simulation 
within the available precision (up to 5%). 

The total efficiency of this selection, defined as the fraction of B° — > j/xpK^ decays produced 
with pj > 5 GeV and |y B | < 2.2 that pass all criteria, is 1.7%. The efficiency ranges from 
0.7% for pj ~ 5 GeV, to 11.4% for pj > 24 GeV, with roughly equal losses due to the dimuon 
kinematic acceptance, the dimuon trigger, and the Kg reconstruction. 

The proper decay length of each selected B° candidate is calculated as ct = (m^/ pj)L xlJ , where 
the transverse decay length L xy is the vector s pointing from the primary vertex [15] to the B° 
vertex projected onto the B° transverse momentum vector: L xy = (s ■ p B ) / | pj | . 

Backgrounds are dominated by prompt and non-prompt J/tp production, with non-prompt con- 
tributions from sources peaking and nonpeaking in m^, as shown in Fig. [I] In particular, mis- 
reconstructed b-hadron decays to final states with a ]/tp, such as B — > J/i/?K*(892), produce 
a broadly peaking structure in the region < 5.2 GeV. A study of the dimuon invariant 
mass distribution confirms that the contamination from events containing a misidentified ]/xp 
is negligible after all selection criteria have been applied. 

The signal yields in each p B and | y B | bin are obtained using an unbinned extended maximum- 
likelihood fit to t«b and ct. The likelihood for event j is obtained by summing the product of 
yield n, and probability density V\ for each of the signal and background hypotheses i. Four 
individual components are considered: signal events, prompt ]/ip events, non-prompt b — > ]/ip 
events that peak in (peaking), and non-prompt b — > ]/tp events that do not peak in 
(nonpeaking). The extended likelihood function is the product of likelihoods for all events: 
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The probability density functions (PDFs), V{, with shape parameters a t - for and /3, for ct, 
are evaluated separately for each of the i fit components. The yields n, are determined by 
maximizing C with respect to the yields and a subset of the PDF parameters. 

The PDF shapes are described below with the parameters obtained from data when possible. 
The niv PDFs are as follows: the sum of two Gaussian functions for the signal; exponential func- 
tions for the prompt and nonpeaking backgrounds; and a sum of three Gaussian functions for 
the peaking background. The resolution on for correctly reconstructed signal events from 
simulation is approximately 20 MeV. The ct PDFs are as follows: a single exponential function 
convolved with the resolution function to describe the signal and peaking background com- 
ponents, where the lifetimes are allowed to be different; the sum of two exponential functions 
convolved with the resolution function for the nonpeaking component; and the pure resolu- 
tion function for the prompt ]/ip component. The resolution function, a sum of two Gaussian 
functions, is common for signal and background and is measured in data to have an average 
resolution of 71 fim. 

The fit proceeds in several steps such that all background shapes are obtained directly from 
data, except for the peaking component which is taken from simulation, as are the signal 
shapes. This technique relies on the assumption that in the region 5.4 < < 5.7 GeV (side- 
band) there are only two contributions: prompt ]/tp and nonpeaking background. To obtain 
the effective lifetime distribution of the nonpeaking background, the and ct distributions 
in the sideband region are fit simultaneously for events in the inclusive B° sample defined 
by pj > 5 GeV and |y B | < 2.2. In the second step, the signal B° lifetime in the inclusive 
sample is determined by fitting ct and simultaneously in the full range. The result, 
cr = 479 ± 22 jim (statistical uncertainty only), is in agreement with the world-average value, 
457 ± 3 Jim fl4). With the effective lifetimes for signal and non-prompt background fixed, the 
signal and background yields are fit in each bin of p B and | y B | , together with the parameters 
describing the ct resolution and the shapes of the prompt and nonpeaking components in m^. 

The accuracy and robustness of the fit strategy were demonstrated by performing a large set of 
pseudo-experiments, with each one corresponding to the yields observed in data, where signal 
and background events were generated randomly from the PDFs in each bin. No significant 
biases were observed on the yields, and the statistical precision of the test was taken as the 
systematic uncertainty due to potential biases in the fit method. The fit uncertainties were also 
observed to be estimated properly. 

The fitted signal yields in each bin of p B and |t/ B | are summarized in Table [lj Figure [l] shows 
the fit projections for and ct from the inclusive sample with p B > 5 GeV and |i/ B | < 2.2. The 
total number of signal events is 809 ± 39, where the uncertainty is statistical only. 

The differential cross section is calculated in bins of p B as 

dajpp -» B°X) = n sig 

&p\ 2-e-B-L-Ap*' K ' 

and similarly for |y B |, where n s j g is the fitted number of signal events in the given bin, e is the 
efficiency for a B° meson to pass all the selection criteria, L is the integrated luminosity, Ap B is 
the bin size, and B is the product of branching fractions B(B° -4 ]/f K^) = (4.36 ± 0.16) x 10~ 4 , 
B(J/ip -> ]i+y--) = (5.93 ± 0.06) x 10~ 2 , and B(K° S n+n-) = 0.6920 ± 0.0005 [14J. The ad- 
ditional factor of two in the denominator accounts for our choice of quoting the cross section 

for B° production only, while n S i g includes both B° and B . The efficiencies are calculated sepa- 
rately for each bin, always considering only mesons produced with |i/ B | < 2.2 (p B > 5 GeV) for 
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Figure 1: Projections of the fit results in (a) and (b) ct for pj > 5 GeV and |y B | < 2.2. The 
curves in each plot are as follows: the sum of all contributions (blue solid line); the prompt ]/ip 
(green dotted); the sum of the prompt ]/ip and peaking background (red dashed), and the sum 
of all backgrounds (purple dot-dashed). 
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Figure 2: Measured differential cross sections (a) da/ dpj and (b) da I dy B compared to the theo- 
retical predictions. The inner error bars correspond to the statistical uncertainties and the outer 
error bars represents the uncorrelated systematic uncertainties added in quadrature to the sta- 
tistical uncertainties. Overall uncertainties of 4% for the luminosity and 3.8% for the branching 
fractions are not shown. The solid and dashed (blue) lines are the MC@NLO prediction and its 
uncertainty respectively. The dotted (red) line is the pythia prediction. 



Pj ( |y B |) bins, and take into account bin-to-bin migrations (< 1%) due to the resolution on the 
measured pj and |y B |. 

The cross section is affected by systematic uncertainties on the signal yield and efficiencies, 
which are uncorrelated bin-to-bin and can affect the shapes of the distributions, and by un- 
certainties on the branching fractions and luminosity, which are common to all bins and only 
affect the overall normalization. The uncertainty on the signal yield arises from potential fit 
biases and imperfect knowledge of the PDF parameters (4-7%), and from effects of final-state 
radiation and mismeasured track momenta on the signal shape in (1%). Uncertainties on 
the efficiencies arise from the trigger (2-3%), muon identification (1%), muon tracking (1%), 
Kg (5%) and B° (3%) candidate selection requirements, acceptance (2-3%), dimuon correlations 
(1-5%) and pj and |i/ B | mismeasurement (1%). The first five efficiency uncertainties are de- 
termined directly from data, while the last three are determined by simulation. The largest 
uncertainties on the efficiency arrise from the Kg reconstruction, which is dominated by the 
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Table 1: Signal yield n S i g , efficiency e, and measured differential cross sections da/dp B and 
dcr/dy B , compared to the MC@NLO [22] and PYTHIA ||16| predictions. The uncertainties in 
the measured cross sections are statistical and systematic, respectively excluding the common 
luminosity (4%) and branching fraction (3.8%) uncertainties. The uncertainties on the signal 
yields are statistical only while those on the efficiencies are systematic. 



P? (GeV) 




e (%) 


da/dp\ (^b/GeV) 


MC@NLO 


PYTHIA 


5- 


10 


240 ± 23 


0.65 ± 0.05 


5.20 ± 0.50 ± 0.59 


3.66 


7.42 


10- 


13 


169 ± 17 


3.32 ± 0.28 


1.196 ±0.121 ±0.117 


1.13 


2.14 


13- 


17 


193 ± 16 


6.37 ±0.51 


0.535 ± 0.045 ± 0.051 


0.49 


0.83 


17- 


24 


138 ± 13 


9.60 ± 0.76 


0.145 ± 0.014 ± 0.014 


0.15 


0.24 


24- 


40 


70 ±9 


11.40 ± 1.04 


0.027 ± 0.003 ± 0.003 


0.025 


0.035 


\v l 


i 


n sig 


e (%) 


da/dy B (fib) 


MC@NLO 


PYTHIA 


0.0- 


0.5 


145 ± 14 


1.34 ±0.10 


7.63 ± 0.74 ± 0.76 


6.21 


12.41 


0.5- 


1.0 


141 ± 15 


1.38 ±0.10 


7.20 ± 0.75 ± 0.71 


6.14 


12.01 


1.0- 


1.4 


167 ±17 


1.93 ±0.15 


7.61 ±0.77 ±0.83 


5.81 


11.24 


1.4- 


1.8 


229 ± 21 


2.51 ±0.21 


8.06 ± 0.74 ± 0.89 


5.38 


10.36 


1.8- 


2.2 


128 ± 17 


1.69 ±0.14 


6.71 ±0.87 ±0.80 


4.81 


9.26 



displaced hadronic track efficiency and is measured by comparing the reconstructed Kg life- 
time with the known value, and the dimuon correlation uncertainty, which is taken as 100% of 
the correction applied to account for the correlations. The difference between the kinematically 
reweighted and unreweighted results (3-5%) is taken as an additional systematic uncertainty. 
The bin-to-bin systematic uncertainty is computed as the sum in quadrature of the individual 
uncertainties, and is summarized in Table [l] In addition, there are normalization uncertainties 
of 4% from the luminosity measurement and of 3.8% from the branching fractions [fl"4|. 

The differential cross sections as functions of p B and |y B | are shown in Fig. [2] and Table [l] 
They are compared to the predictions of MC@NLO ||22| using a b-quark mass of 4.75 GeV, 

renormalization and factorization scales fi = \Jm^ + p\, and the CTEQ6M parton distribution 
functions fl23| . The uncertainty on the predicted cross section is calculated by independently 
varying the renormalization and factorization scales by factors of two, by ±0.25 GeV, and by 
using the CTEQ6.6 parton distribution functions. For reference, the prediction of PYTHIA [16] is 
also included, using a b-quark mass of 4.80 GeV, CTEQ6L1 parton distribution functions p3| , 
and the Z2 tune [24] to simulate the underlying event. The measured pi spectrum falls slightly 
faster than predicted by MC@NLO, while the y spectrum is measured to be flatter than the 
PYTHIA prediction and in agreement with the MC@NLO prediction within uncertainties. The 
integrated cross section for p B > 5 GeV and \y B | < 2.2 is calculated as the sum over all p T bins, 
without an upper limit for the highest pj bin, to be 33.2 ± 2.5 ± 3.5 fib, where the first uncer- 
tainty is statistical and the second is systematic. The result is compatible with the prediction 
from MC@NLO (25.2+g| ar >d below the prediction from PYTHIA (49.1 fib). 

In summary, the first measurements of the differential cross sections da/ dp B and da I dy B for 
mesons produced in pp collisions at ^fs = 7TeV have been presented using the decay 
B° — > J/]/? Kg. The measurements cover a range in p B from 5 GeV to more than 30 GeV, and 
the rapidity range \y B \ < 2.2. The total cross section in this kinematic region lies between 
the central values of the MC@NLO and PYTHIA predictions, with a rapidity distribution that is 
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flatter than PYTHIA. It is also in agreement within uncertainties with the measured B + cross 
section 0. 
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